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Abstract 

A non-equilibrium plasma was studied using classical electrodynamic field theory. Non-linear 
interaction terms contribute to a finite lifetime for the dressed electrodynamic field. The lifetime 
exhibits a ~ n~^Te^'^T~'^Tr^'^ dependence, where n is the number density, Tg is the electron tem- 
perature, Ti is the ion temperature, and T^. is the temperature of the radiation field. The resulting 
width of the plasmon resonance is shown to decrease as equilibrium is approached. Dynamic 
screening leads to opaqueness of the plasma for low energy electromagnetic radiation. This leads 
to a quadratic correction to the quartic Stefan-Boltzmann law. We also briefly discuss the effect 
of dynamic screening on fusion rates. Solitonic solutions to our non-linear wave equation allow 
localization of positive charges, which may enhance fusion rates. 
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I. INTRODUCTION 



Non-equilibrium plasmas have been investigated theoretically in several regimes by many 
authorsii^iM. These papers focus on approximations to deal with inter-particle collisions 
in plasmas. They develop various models to describe the temporal progress of a non-local- 
thermal-equilibrium (NLTE) plasma towards equilibrium. 

Recently, probative experiments have ben performed on the passage towards equilibrium 
of a plasma in non-local thermal equilibrium (NLTE)^. The paper by Taccetti at al^ gives a 
putative bound of about Ins as the time required for the NLTE plasma in their experiment 
to achieve equilibrium. There will be a similar bound for any other NLTE plasma. On 
time scales much smaller than this measure, one may safely assume that a NLTE plasma 
is frozen, viz., the electron temperature is different from the ion temperature. But light 
can travel appreciable distances on this time scale (e.g. 3cm in 0.1ns), so that we need to 
consider the time dependence of electromagnetic fields. By extending the static methods of 
Jaffe and Brown^ and Chitanvis^, we can consider the classical statistics of the electrody- 
namic field in a non-relativistic NLTE plasma. By varying the electron temperature, the 
ion temperature, as well as the temperature of the electromagnetic field, we can study para- 
metrically the passage of an NLTE towards local thermal equilibrium (LTE). Phenomena 
such as Bremsstrahlung and inverse Compton scattering of hot electrons which equilibrate 
the temperature of electrons and the radiation field, and electron-ion collisions, which equi- 
librate the energy between electrons and ions could be considered in the future to describe 
the temporal development of a non-equilibrium plasma towards LTE. 

We have focused here on the dynamic dielectric constant of a non-equilibrium plasma, 
which can be measured in a laboratory. The dispersive properties of this quantity will be 
discussed. We propose that experiments be performed to test our predictions. 

II. FORMULATION OF THE STATISTICAL THEORY 

The electrostatic interaction between any two charges sitting in a collection of other 
charges is inevitably screened by intervening charges. This effect becomes more dominant 
as the density of charges increases. Thus calculating the screening of the bare Coulomb 
interaction could be important for high-density plasmas. The paper by Varela et al^ gives a 



recent overview of various methods brought to bear on this aspect of plasma physics. Some 
of these methods develop ideas based on the work of Fisher and Levin^. Our field-theoretic 
method bears a close resemblance to the approach of Brown and Jaffe^ and Chitanvis''. We 
shall focus on the calculation of non-linear effects on the electrodynamics of a NLTE plasma. 
These non-linear effects are expected to dominate at high charge densities. 

An assumption we will make is that on the scale that fluid motion takes place, the 
individual plasma components have separately reached thermal equilibrium, such that the 
ions of a specie possess a temperature which is different than the electron temperature. 

The scalar potential of a plasma composed of a single atomic specie, having an atomic 
number Z, and number density n is governed by: 



c-'dt<p - V> = 47rp 

P = P++ P- 
p+ = e n Z exp{— Zecj)/ ksTi) 

p_ = —e n Z exp{e(f) / kTf.) (1) 

where c is the speed of light, p is the charge density, where ks is the Boltzmann constant, e 
is the electronic charge, Tj is the ion temperature, Tg is the electron temperature, and n is 
the average number density of the plasma. Since we consider a non-relativistic plasma, we 
can neglect the effects of the vector potential, which are 0{v/c), compared to the potential 
(f). We have assumed that positive and negative charge densities are governed by Boltzmann 
distributions, and the electron temperature is in general different from the ion temperature. 
In so doing, we are making the assumption that quantum effects are unimportant. Now, for 
the case of plasmas in local thermal equilibrium, Gruzinov and Bahcall^, among others, have 
estimated that quantum effects give rise to small deviations from Eqn^ for the screening 
length. This could be important for certain applications. 

Expanding the right hand side to third order in the electrostatic potential (basically in 
powers of a parameter F which is analogous to the usual plasma expansion parameter), it 
can be shown that: 

c-^aff - « (^^^) (-i + "ii' - i^') (2) 
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where r = Te/Ti, F is a parameter analogous to the usual plasma parameter. Basically, 
we use A as a length scale rather than the mean free distance, since A appears naturally 
in the theory, rio is the average number density of the plasma, and Xdh is the standard 
Debye-Huckel screening length, which holds in the case of a single-component plasma. 

Equation El demonstrates as long as T « 1 the coupling constants a2 < I, < 1, as 
long as Zt not too high. In this sense, we have tried to ensure that our series expansion is 
at least formally convergent. The value of r ranges between a number much greater than 
one to unity, as the non-equilibrium plasma tends towards equilibrium. As an example, 
we note that at Tg = IkeV, A ~ 10~^cm, F ~ 0.06. Hence, our expansion is valid for 
Zt < 1/0.06 ~ 16.7. For Z ~ 1, our expansion is applicable to plasmas fairly out of 
equilibrium. 

Furthermore, note that retaining non-linear terms on the right hand side of Eqn^ in 
powers of F, is equivalent to seeking corrections due to the high density of the plasma^. The 
number of terms that need to be retained thus depend on F. 

Now in the linear case when o"2 = o"3 = in EqnEl it is easy to see that our plasma has an 
effective screening length A which is given by the standard single- component Debye-Huckel 



length divided by a factor of a/ {Zt + 1)Z. Thus this correction can be significant for large 
Z, and for cases when there is a large temperature difference between the electrons and ions. 

Note that EqnlHJfor the electrostatic potential can be obtained by extremizing the fol- 
lowing Lagrangian density: 



m =\{{dt4Y-\v'^\')-v{4>) 



(4) 



wher t' = ct/A, where the prime on the gradient operator indicates we are using A as a 
length scale. 

The cubic term can be interpreted as a three- wave interaction, and the quartic a four- wave 
process- If second quantization is employed, these terms would correspond to three-photon 
and four-photon processes. One may perhaps expect these terms to be important in the 
high energy density regime. Moreover, the theory bears a formal resemblance to the theory 
of phase transitions. However, we do not have a phase transition to contend with, and hence 
the effect of fluctuations could be minimal, and mean field theory (EqnlHl corresponds to a 
dynamic Debye-Huckel approximation) would be sufficient to describe the system. In fact 
detailed calculations presented in this paper show that higher order processes are not be 
important as far as renormalization of the plasmon frequency are concerned. But they are 
essential if one wishes to consider the lifetime of the dressed electromagnetic field. 

Since we have a time-independent Lagrangian density, the momentum vr = dC/d(j) = dt'ip- 
Hence the Hamiltonian density ( i.e. the energy density) is: 



Note that the Hamiltonian density (energy) is positive semi-definite. We define the following 
partition function in order to study the statistical mechanics of the system described by this 
energy density: 



where T^/Tg is the temperature of the electromagnetic field in units of the electron temper- 
ature, since we are using ksTf. as the unit of energy. The functional integral over all fields 
denotes a sum over all thermodynamic states. In this sense, we take into account fluctua- 
tions around Eqn. IHl We consider Eqn. 121 a mean-field approximation, and a non-linear, 
dynamic extension of the standard Debye-Huckel equation. 
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In order to simplify calculations, we will define $ = {^/TJT\.) 0, and 02 ^ ci'2 = 
^/TrjTf. (72, and o"3 — CTg = {Tr/Te) (73. Sincc we generally expect Tg >> for an NLTE 
plasma, ag << a^. 

In terms of $, the partition function can be written as: 



Z = j V(^exp{- j d^x' ni^)) (7) 

We can now use standard techniques from (Euclidean) field theory^^ to obtain the lowest 
order non-trivial corrections to the self-energy, or equivalently, the screening length. The 
correction to 0{a2) is obtained by noting that when higher order corrections due to the 
non-linear term in V^($) is taken into account: 

Goiuj,k) ^ (G,\uj,k) ~tiu;,k)y' (8) 

where Go{uj, k) denotes the momentum-space inverse of the operator {—d'ji — + 1) in 
the units indicated above, viz., GQ{uj,k) = (cu^ + k"^ + and where f]{uj,k) is the self- 
energy. Note that this is somewhat different than the usual Klein-Gordon-like operator in 
momentum-space, (a;^ — k"^ — 1)^^. One can obtain an equivalence by letting k ^ ik and 
m im. 

We used Mathematica to generate symbolically the second order contribution (the first 
non-trivial correction) from the cubic term in the energy functional. The net result is: 

'^polarization-likei^iP) = ^ j ^^^^4 ^o(^7 k)GQ{uJ — Vl,p — k) 



Go{u;,p) = (9) 



The form of the integral can be evaluated using dimensional regularization^i: 
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48-167r2^ V / V /V y^l-4/(Lu2 + P) - 1^ 

We will work in the regime where Te >> T^, so that the quartic term in \^($) can be 
safely ignored, using the scaling arguments described earlier. 



The appearance of an imaginary term preceding the logarithm gives rise to a hfetime 
for electromagnetic excitations, analogous to the lifetime of a quasi-particle. Furthermore, 
there is an additional contribution to the lifetime for 1 < 4/(c(j^ + k"^). 

The net effect is that we account to 0{a2) the screening effects on the pure Coulomb 
potential. It is straightforward to obtain the effective wave-equation, obtained in lieu of 
EqnH 



The effective wave equation follows simply, by extremizing this Lagrangian, and the 
spectral form of its Green's function can be obtained as follows: 

Q{UJ, k) = (u'^ - k'^ -1 + f^polarization-HkeiuJ, k))^ (12) 

Basically, Eqn. yields a plasmon resonance at : 

uj = ±^c^k^ + ujl (13) 

where c = 1 in the units we are using. This is a perturbative result, valid for o"2 < 1- O^i^ 
might be tempted to assume from Eqns. ^1 that the overall contribution is << 1, given 
the small numerical pre-factor. As speculated earlier, this is indeed the far as the 

plasmon frequency is concerned- In the units chosen, A = 1, and tu^ ~ c^/A^ + 0{(J2 ) ~ 1. 
The effects of interaction on the renormalization of the plasmon frequency can be gauged 
by studying the zero of Re{Q^^{uj, k = 0.1)) as displayed in Fig. [H The zero occurs at 
the plasmon frequency Up. The plasmon frequency is insensitive to the coupling constants, 
which in turn depend on the parameters of the plasma. 

The negative of the imaginary part of the structure factor Q{uj,k) yields the spectral 
weight. The spectral weight manifests the plasmon resonance. This quantity should be 
accessible to experiments. Qualitatively, this is what one expects to occur in a plasmap^. 
For sufficiently large frequencies, far from the plasmon frequency, one expects a long-lived 
wave. This effect is exemplified in Fig|21 This figure shows the spectral weight for two 
different values of the coupling constant. The red curve refers to = 1 and signifies a 
case more out of equilibrium, and the blue curve with = 1/2 refers to case closer to 
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FIG. 1: The blue curve corresponds to CJ2 = 1 and has a zero at the plasmon frequency of « 1 
(in the units chosen), while the red curve corresponds to a'2 = 10 and also has a zero at ~ 1. 
This shows that the plasmon frequency does not depend sensitively on the coupling constants, 
which in turn depend on the plasma parameters. Note that while the theory itself is formally 
valid for values of the coupling constant less than one, we have used larger values of the coupling 
constants only for exposition. 

equilibrium. The spectral weight is sensitive to the value of the coupling constants (see 
Figure Ej) over part of the frequency spectrum. Notice also that the width narrows as the 
system approached equilibrium. 

III. DYNAMIC SCREENING EFFECTS 

The dielectric constant e{uj,k) of the plasma may be obtained via the following formal 
identification: 
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FIG. 2: This plot of the spectral weight. It clearly shows the plasmon structure, broadened by 
three-wave interactions. The blue curve corresponds to (T'2 = 1/2 and refers to a case closer to 
equilibrium, and the red curve corresponds to an NLTE case of = 1. 



Q{uJ,k) — yuj'^ — /C^ — 1 + T,polarization-like{^ , k))j 
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e{uj, k) =1 ^ (14) 

We have plotted in Fig. Elthe negative of the imaginary part of the inverse of the dielectric 
constant, which is proportional to the structure factor S{uj, k)^- of the NLTE plasma. The 
structure factor may be probed in future experiments. 

IV. RELAXATION RATES 

Third-order interactions (corresponding to three photon processes) lead to finite lifetimes 
for electromagnetic excitations. Jaffe and Brown^ have pointed out in the static case, there is 
a correspondence between our approach, which focuses on the electromagnetic field, and the 
ensemble consisting of electrons and ions. A similar connection must exist for the dynamic 
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FIG. 3: This plot of the negative of the imaginary part of the inverse of the dielectric constant, 
which is proportional to the structure factor. It clearly shows the plasmon structure, broadened 
by three- wave interactions. The blue curve corresponds to = 1/2 and refers to a case closer to 
equilibrium, and the red curve corresponds to an NLTE case of = 1. 

case considered here as well. Relaxation rates for collisions between ions and electrons can 
be computed for various regimesi^. Not only does our theory account for such inter-particle 
collisions in an integrated sensed, but addresses radiative processes as well. The inverse of 
the rate we consider here is a measure of the time required for electromagnetic energy in a 
specific part of the spectrum to be given to other parts of the spectrum. 

To the first order in the coupling constant cTg, the relaxation rate, written as an inverse 
lifetime is: 




^™ - (x) [ii) 



for 1 > 4/(^17^ + fc^). Here we have re-introduced dimensionful parameterization. 

For the case of extreme non-equilibrium, for Z >> 1, when r >> 1, we can extract the 



main functional dependencies as: 
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This result can be compared to the usual Spitzer rate (two-temperature plasma), which 
is restricted to electron-ion collisions, and yields a Te temperature dependence. 

For 1 < 4/(a;^ + ^^), we obtain an additional contribution to the lifetime which depends 
on the region of the spectrum being probed: 
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V. RADIATIVE PROPERTIES OF A NLTE PLASMA 

We have found that the dispersion relation is given approximately oj = ^/c^ + Up. This 
implies that electromagnetic waves not satisfying this relation exist at best in an evanescent 
fashion in the plasma. One expects on physical grounds that this absorptive-transparency 
phenomenon will affect the radiative properties of a black-body in which a plasma exists as 
well. 

One can estimate this effect by modifying the standard calculation of the energy density 
of the electromagnetic field by noting that the density of states g{u!) per unit volume for the 

present case is: 



g{iS)dijj — dcu Vo; > cUp (18) 

Here we have restored dimensionality of all parameters. The resulting energy radiation 
fiowing per unit area, per unit time, to leading order in cUp is then: 

S{Tr,u;p) = asBT^ ( 1 - ^^^^ + ... ) (19) 



where gsb is the Stefan-Boltzmann constant. 
We note the following estimates: 

• The plasmon frequency uj^ ~ IQ^^Hz at an electron temperature of IkeV, and a number 

density ^ 10^^cm~^. This corresponds to photon energy huj^ ~ lOey. This suggests that our 
corrections to the Stefan-Boltzmann law can be significant for a radiation field temperature 
Tr - lOeV. 



VI. QUANTUM TUNNELING 

We will consider briefly in this section the effect of dynamic screening on fusion rates. 
Toward the end of the section, we indulge in speculating how the thermal rate of fusion 
may be aff'ected in NLTE plasmas. This may be considered a prelude to a full-fledged study 
of nuclear fusion rates in non-equilibrium plasmas. We have restored dimensionality of all 
parameters in this section. 

Now the WKB penetration factor V{E) for charged nuclear particles interacting via a 
pure Coulomb term at sufficiently large distances compared to the nuclear lengths is: 



V{E) = exp (^-2 K{r)dr^ 
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where ri — R is the turning point where nuclear attractive forces begin to dominate 
Coulomb repulsion, and the other turning point r2 = ZZ'e^/E, Ze and Z'e being the 
charges on the interacting particles, n is the reduced mass of the two nuclear interacting 
particles, v is the relative speed, and E is the sum of the kinetic energy and the static 
screening effect. 



Now, from Eqn IlUl the pure Coulomb interaction is modified thusly: 



1 exp(-r Jujl/c^-El/{c^h')) 

- ^ (21) 

where for sufficiently low values of the kinetic energy E^, the radical in Eqn|^ remains real. 
It is seen that accounting for the finite speed of light causes the screening length to increase, 
i.e., the range of the interaction increases, allowing for repulsion over a larger distance. 
Hence we expect the penetration factor to decrease. 

Upon expanding in powers of r, since we are only interested in short separations, and 
further, retaining only leading order terms in E^, we obtain an addition to the argument in 
the penetration factor: 
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a = ZZ'e' (22) 

Hence, by accounting for the finite speed of light, we have shown that retardation effects 
will affect the fusion rate. Salpeter— referred to dynamic effects in passing in his seminal 
paper, while more recently Shaviv and Shavi\»i^, and Bahcall et a^^ have considered such 
effects in more detail. Here we have presented an alternative treatment of the same effect. 

Upon utilizing Up/c ~ lO^cm"^, Ek ~ lOkeV, R ~ lO^^^cm, /i ~ 10~^^g, we find that the 
first term in Eqn|221is greater than the first, so that the net effect for this set of parameters 
is to decrease the fusion rate, albeit by a small amount 0(10"'^), compared to Eqn l2Ul which 
yields a contribution of ~ C(10^^). 

Note that this theory is based on the assumption that the charge density follows the 
classical Boltzmann distribution even at nuclear distances. This point has been questioned 
by many authors, e.g. Gruzinov and Bahcall^^ for the case of a plasma in local thermal 
equilibrium. They find that quantum effects change the screening length by approximately 
1% in the solar interior. We will consider such effects for a NLTE plasma in a future paper. 

It is well-known that the fusion rate at a finite temperature is directly proportional to 
the relative fraction of the reactants. Hence it is possible that if concentration of reactants 
is increased locally, it would enhance fusion in that spatial region. To see that such charge 



configurations may occur in our system, let us consider the wave equation to the second 
order in the potential: 



(23) 



We seek a traveling wave solution in one dimension of the form $(x — ut), where the 
wave-speed u is measured in units of the speed of light. Redefining ^ = {x — ut)/yl — v?, 
with M < 1, we see that for the special case of 2; = 0, the solution is: 



This is of course a traveling soliton. It is interesting to see what charge distribution 
is carried along with this soliton. This can be estimated by computing the curvature of 
the solution, viz., — (9|<l>(.^). We illustrate this is FigEJ Notice that the blue curve is the 
potential, which shows a change in the curvature as we move away from the origin. This is 
manifested in the charge density, displayed as the red curve. Notice how the positive charge 
gets localized by a surrounding negative charge density. Furthermore, we have lowered in 
Fig. I^he charge density by a scaling factor oi a'2 <1. This represents a deviation from the 
average charge density, which will be much lower. The fact that positive charges (nuclei) 
have been localized will enhance the rate of fusion. Of course, whether such localization 
occurs in a given plasma depends on initial and boundary conditions on the NLTE plasma. 

VII. CONCLUSIONS 

We studied parametrically dynamic screening effects in a three-temperature plasma. A 
statistical theory of the electromagnetic field interacting with charges was developed for 
this purpose. Diagrammatic techniques were utilized to obtain leading order estimates of 
the lifetime of the plasmon resonance due to three-wave interactions. Effects of dynamic 
screening on radiative properties were deduced. We speculated on the manner in which 
solitonic solutions could allow localized enhancement of positive charges, which in turn 
would increase fusion rates. 

We propose that the techniques developed, and results obtained in this paper could be 
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FIG. 4: The blue curve is a solitonic solution to Eqn l23l The red curve is the corresponding charge 
density. Note how non-linear effects serve to localize the positive charge. 

tested in future experiments. 
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